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Experiments strongly suggest that the flavour mixing responsible for the atmospheric 
\ neutrino anomaly is very close to being maximal. Thus, it is of great theoretical as well 

Q\ ' as experimental importance to measure any possible deviation from maximality. In this 

context, we reexamine the effects of matter interactions in long baseline neutrino oscilla- 
i , tion experiments. Contrary to popular belief, the muon neutrino survival probability is 

J> \ shown to be quite sensitive to matter effects. Moreover, for moderately long baselines, 

the difference between the survival probilities for and is shown to be large and 
. sensitive to the deviation of {U^l from maximality. Performing a realistic analysis, we 

\ demonstrate that a muon-storage ring i/-source alongwith an iron calorimeter detector 

can measure such deviations. (Contrary to recent claims, this is not so for the NuMI- 
^ > ■ minos experiment.) We also discuss the possible correlation in measuring and U e z 

in such experiment. 

a: 

1 Introduction 

Neutrino masses and mixings continue to intrigue us, despite the continued efforts of many 
recent experiments [1-6]. We do know that the muon- neutrino (z/ M ) mixes almost maximally 
(sin 2 29^ 3 > 0.92 at 90% C.L) with another species z/ 3 (which could be identified with the tau- 
neutrino, v T ) and that the mass splitting |5m 2 2 | ~ 2 x 10~ 3 eV 2 . Furthermore, the electron- 
neutrino {y e ) mixes with a combination of and v T with a similarly large angle 9 e 2 — 30° 
but a far smaller mass splitting (5m|i ~ 7 x 10 -5 eV 2 ). And, finally, the third mixing angle 
#i3 is constrained, by the chooz experiment [3], to be small. What we do not know though 
are the sign of <5m 2 2 , the precise value of 0^ and the extent to which (9 m3 may, if at all, differ 
from maximality. Since such knowledge is obviously essential to the formulation of a theory of 
neutrino flavours, the importance of precise determinations cannot but be underemphasised. 
Neutrino oscillations provide just the proper platform for such an exercise. 

As is well known, during passage through dense matter, various amplitudes for neutrino 
oscillations may be magnified to a great degree [7] thereby rendering them measurable with 
relative ease. While propagating in a medium, all the three neutrino flavours interact with 
matter via neutral curent (NC) interactions with equal strength. Those of the electron flavour, 

pi , have, in addition, a extra charged current interaction as well. Since it is this extra 
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interaction that provides flavour-dependence, all matter-induced effects in neutrino oscillations 
(including those in the v^-v T sector) are proportional to the size of the relevant mixing with 
v e . 

Recent theoretical advances in the quest of producing very high intensity muon sources 
engender optimism regarding a future neutrino factory wherein an intense beam of muons 
is to be accelerated to a not too high energy and stored in a storage ring with a straight 
section directed towards a neutrino detector. Muons would decay in this straight section 
thereby producing high intensity neutrino beams (both of electron- and muon-types) that are 
highly collimated in the direction of the decaying muons [8]. The efficacy of muon storage ring 
neutrino sources in performing precision measurements of neutrino parameters has already been 
discussed in the literature [9]. The advantages over neutrino experiments with conventional 
neutrino beams (arising from 7r ± decay) are manifold: (i) a precise knowledge of the and 
v e fluxes helps reduce the systematic errors; (ii) assuming a /x-beam of, say, 20 GeV energy, 
such neutrinos have, on the average, energies higher than those of conventional neutrino beams 
thereby increasing the v cross-section at the detector; and (in) the aforementioned collimation. 

Previous studies of the matter effect in Ref. [10] at a long base line experiment were con- 
centrated mainly on the measurement of the rate of wrong sign muon events, which, in turn is 
proportional to the transition probablity P Ve ^ v . Such a measurement would enable a deter- 
mination of (or constraining of) sin#i 3 provided independent and precise measurements of the 
mixing angle 6*23 and mass-square difference Sm^ exist. Possible measurements of neutrino 
mixing parameters from long baseline experiments using conventional neutrino beams as well 
as from a neutrino factory, and a water Cerenkov detector has also been discussed in Ref. [11]. 

In the present article, we demonstrate that a sizable matter effect is also evinced by a 
measurement of the survival probability. In addition, such a measurement has the poten- 
tial of detecting possible deviations from maximality of the mixing angle responsible for the 
atmospheric neutrino anomaly. In a recent study [12], the presence of large matter effect in 

— v T oscillation as well as in survival probability has been stressed. Although there is a 
large matter effect present in — v T sector, measurement of the same is difficult due to both 
the lower detection efficiency for the tau-events as well as the need for a more sophisticated 
detector. 

Whereas much of our analysis would be source and detector independent, we also choose to 
examine the feasibility of making such meaurements in the context of a realistic experimental 
setup. To be concrete, we choose a 50 kT Iron detector, with detection and charge discrimina- 
tion capability for muons provided by a magnetic field. Such a detector was proposed for Gran 
Sasso (monolith) [13] and, more recently, for a location in India (INO) [14]. The latter is 
contemplated primarily as a detector for atmospheric neutrinos and also as an end detector for 
a future neutrino factory beam. A muon detection threshold of 2 GeV [14] has been used in our 
calculation. As for the baseline, we illustrate our results for two particular choices. The first 
one coresponds to a distance from JHF to INO (~ 5000 Km) and the other one corresponds to 
even a longer baseline, from FERMILAB to INO (~ 10000 Km). We will also compare results 
with a baseline of 732 Km, corresponding to the Fermilab-MINOS distance. 

The plan for the rest of the article is as follows. In section 2, we discuss neutrino oscillation 
in the presence of matter. We also define the asymmetry of /x + and fi~ events and explain its 
correlation to U e ^ and U^. We present the expressions for the asymmetry for two different 

aproximations: when matter effect is small compared to i.e for short base line experiments 
like MINOS [15] and also when the matter term is comparable to the i.e., for long base line 
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experiments. In section 3, we briefly review and reexamine the proposal, in a recent article [15], 
to measure matter effect and the above correlations in the context of the MINOS experiment. 
Section 4 deals with the expectations at lomg baseline experiments, both for the survival 
probability as well as for the v e to transition probability. An analysis of the findings in the 
context of a realistic experimental situation is effected in Section 5. Finally, we conclude in 
Section 6. 



2 Neutrino Oscillation in presence of matter 

We begin by reviewing the passage of neutrino through matter. While analytical results are 
possible, the exact answers are cumbersome and not very useful. However, it is useful to 
consider approximations such as a constant density profile for the earth. While such an as- 
sumption is obviously not a very accurate one, especially when neutrinos traverse a longer 
distance through the earth, the approximation serves to give us some physical insights to the 
problem we are dealing with. Furthermore, results such as those we shall obtain are still appli- 
cable if the matter density within the earth may be approximated to be piecewise constant. It 
must be kept in mind though that the quantitative results presented in this article are obtained 
not by means of such an approximation, but by explicitly accounting for the varying density 
profile of the earth using the Preliminary Reference Earth Model (PREM) [16]. 

As with the quark sector, neutrino flavour states can be expressed in terms of mass eigen- 
states through a relation of the form 



W, 



e) = J2 U ™ \ V i) ' W 



where U is a 3 x 3 unitary matrix known as the Pontecorvo-Maki-Nakag awa-Sakata mixing ma- 
trix. In the absence of Majorana phases (assumed to be so henceforth) 3 , U can be parametrized 
in terms of three mixing angles and a CP- violating phase <5cp, viz. 



( c 12 c 13 s 12 ci 3 si3e" i<5cp \ 

U = -C 23 Sl2 - S23Sl3Cl 2 e l<5cp C23C12 - S23Sl3Sl2e l<5cp S23C13 

\ S23S12 - C23Sl3Cl 2 e l<5cp -S23C12 - C 23 Sl3Sl2e l<5cp C23C13 / 



(2) 



with Cij = cos 9ij and s^- = sin%. In vacuum, neutrino oscillations are then governed by U 
and the two independent differences in the squares of the neutrino masses, Smfj = mf — to 2 . 
In the rest of our analysis, we will make a further simplifying assumption of ^cp — 0. With the 
above parametrization of U, experimental results can be used to constrain the mixing angles 
9ij and the mass-square differences. Measurements of the solar neutrino fluxes, augmented by 
the KamLAND reactor data give 

0.30 < tan 2 12 < 0.64, 5.4 x 10~ 5 eV 2 < S 21 < 9.4 x 10~ 5 eV 2 (3) 

with the best-fit values given by tan 2 12 = 0.40 and Sm^ = 6.9 x 10~ 5 eV 2 . Experiments with 
atmospheric neutrinos (as well as the K2K experiment) suggest, on the other hand, that 

sin 2 2^23 > 0.86, 1.4 x 10~ 3 eV 2 < |5m 2 2 | < 5.1 x 10~ 3 eV 2 , (4) 



3 Any Majorana phase in the PMNS matrix cannot be probed by an oscillation experiment. 
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with the best-fit value being \8m\ 2 \ = 2.0 x 10 3 eV 2 . And, finally, the negative result of the 
CHOOZ experiment implies that 

sin 2 13 < 5 x 10~ 2 at 99.73% C.L. (5) 

Note that, 

• unlike in the case for Sm^, the sign of 5m| 2 is undetermined; 

• since |5m| 2 | ~ l^iil \^ m 2i\ 4 > the last-mentioned (as well as the solar mixing angle 
612) plays only a subservient role in oscillations; 

• while is constrained to be close to l/y/2, it is still allowed to be non- maximal. 

Working in the limit of a vanishingly small Sm^, all the expressions for the oscilla- 
tion/survival probabilities can be expressed in terms of just two of the mixing angles 5 , viz. 
623 and 0i3 . For a (P M ) of energy E v travelling a distance L in vacuum, the survival proba- 
bility is then given by 

P£ = P m = 1 - 4 IM 2 (1 " IM 2 ) sin 2 (AsiL/2) , (6) 

where, 

U e3 = S13 , U^ 3 = C13 S 2 3 , A31 = — . (7) 

In the presence of matter, the electron (anti-)neutrino sees an additional potential correspond- 
ing to charged current interactions, and the effective Hamiltonian is given by 



(8) 



where A = ±\^2GfN £ for j/J , with iV e being the instantaneous electron density of the matter. 

The survival probablities suffer a consequent modification. While the full expressions are 
cumbersome, it is instructive to look at a simplified form obtainable for a matter density small 
enough that A can be treated as a perturbation in the effective Hamiltonian. In this limit, 
following Ref. [15], 

P^m * ± (AP^/2) + 0(A 2 ) 

(9) 

AP W ~ 4A|f/ e3 | 2 |^3| 2 (l-2|^ 3 | 2 ) 4sin 2 ^-^sin(2tf) + 0(A 2 ) 









0\ 


Ms 













0; 



where 



A^A and * = (10) 

^31 * 



4 From now on we will denote this large mass-square difference by Sm^i 

5 While the analytical expressions below have been derived under this approximation, all of the numerical 
analysis has been performed using the full expression for U and a non-zero value for 612- 
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The asymmetry AP^ is a measure of the matter effect and that it should be proportional 
to both A and |?7 M 3| 2 is obvious. Note also that the has the same interaction with matter 
as the v T with which it primarily mixes and any matter effect can only seep in through a 
mixing with the v e . This, then, accounts for the overall factor of |£/ e 3| 2 . All of these three 
proportionalities quite independent of the approximation, and, in fact, are exact results. The 
factor (f — 2|£7 M 3| 2 ), however, is but a consequence of the approximation of a small matter term 
A and is applicable only for neutrinos traversing small base lines. Since this factor vanishes 
identically for a maximal mixing in the — v T sector, this matter-induced asymmetry is 
potentially a sensitive probe of the maximality of U^. This has already been pointed out [15] 
and examined in the context of the MINOS experiment. In the next section we will briefly 
reexamine their claims. 

Convoluting eqn.(9) with the independent flux and the detection efficiency of a given 
detector, one would obtain the number of events and thereby the event asymmetry 



AN 



N^-N^ = N_-N, 



(11) 



The latter could then be used to place contraints in the | j — j C/" e 3 ] plane. Of course, in a 
realistic case, the effects of a varying matter density has to be taken into account as well and 
this we do include in our analysis. With the effect being proportional to A, it will be magnified 
when the (anti-) neutrino beam traverses regions of high density, or in other words, the core of 
the earth. Note, however, that a large A results in a breakdown of the simplified form given in 
eq.(9) and that the full expression (or, at least, a different approximation) needs to be used. 

The exact expression (for an arbitrarily large but constant density) for the survival proba- 
bility is given by 



1-44 



2 2 -2 
C 23 S 9 m sm 



, z z 



A 31 L 



(l+A-V 



sm 



A 31 LV 



+ c 23 c em sin 



A 3l L 



(l + A + v) 



(12) 



where 



V = + A 2 -2Acos29 13 , and 0„ 



1 

- tan 
2 



-1 



sin 2#i3 
cos 26»i3 - A 



(13) 



We now perform a double expansion 6 in \U e3 \ as well as (3 = \ — while allowing for 

any value for A. Since both the expansion parameters are small, we may retain terms only 
upto, say, the fourth order without any loss of accuracy. Then, for the survival probability, we 



3 This is similar in spirit to, though not the same as, the expansion performed in Ref. [17]. 
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have 



= (l -4/3 2 ) 4 



+ 



+ 



\U e s\ 2 



\Ue3\ 4 

(I -A) 



|4-4-4-^ (i-^) *s 2 *}-4/?{(i-a) 2 4-4 

+4/3 2 1(1 - AA+ 2A 2 ) 4 + 4 - 4 + | (l - A) V S2 * j 

(l-6A + l 4 -4l 3 + 5l 2 ) 4 + (2l + i 2 )4 
-ft 2 (1 - A) 2 c 2 * + (-1 + 6 A - A 2 ) 4 



+ (l - A 2 ) n s 2 * + (l - A) Q (s 2n + 2 s 2e ) 
+ 0(\U e3 \ 4 f3,\Ue3\ 3 (3 2 ,\U e3 \ 2 (3 3 ,\U e3 \(3 A ) , 
where s a = sin(a), c a = cos(o;) with 

A 31 LA AL 



n = 



A * and = ft - * . 



With the survival probability for the antineutrinos being obtained by the replacement A 
in eq.(14), we have, for the probability asymmetry 



-(l-A 2 ) 3 AP 
4l|£/ e3 | 2 |M 2 



(1 - A 2 ) 



444 + (i-^)i% - 



l + A 2 



2A 



+ 2f3{l-A 2 ) 



l + A 2 



+ \U, 



e3 | 



- 4 4 4 + (1 - A ) * s 2 * + ~ s 2 * s 2 ^ 

^ ./i 

T A |i {(l + llI 2 + 2l 4 )4c 2 +l 2 (3-I 2 )44} 
+^ {2 A 2 (1 - I 2 ) 2 $c 2 $-(l + 6l 2 + I 4 ) c*} 



+s 2 * |2(l + 3A 2 )^c 2n - 

-(3 + I 2 )fts 2f , (3-4 4) 
+ 0(|C7 e3 | 2 /3, |t/ e3 |/3 2 ,/? 3 ) . 



1 - 16 A 2 - A 4 



4A 



■ S 2 Q 



Note that 



(14) 



(15) 
• -A 



(16) 



• The asymmetry AP continues to be proportional to each of A, \U e3 \ 2 and {U^ 2 , but not 
to(l-2|^ 3 | 2 ). 

• lim^ s^i/A = \1/ and hence the right hand side of eq.(16) is finite in this limit. 

• To 0(A), the expressions for both P m and AP coincide with those of eq.(9). 
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• The strength of the (9(1) term on the r.h.s. of eq.(16) — which vanishes as A 2 for small 
A — gives a measure of the violation of the (1 — 2|L^ 3 | 2 ) proportionality of AP. 

• While we have performed the expansion upto C ( | t/ e 3 1 4 ) , for all practical purposes it 
suffices to consider only upto 0(\U e3 \ 2 ). 

• The apparent singularities in eqs.(14 & 16) are not real ones with 
limP w = c| + 4/3 2 4 



A->1 



+ \u e3 

+ lt/ e3 l 4 * 2 



(s 2 *-2*4)+4/3 (^ 2 -4) +4/3 2 {2 (l-^ 2 ) sl-Vs^} 



\1> 2 4 

— (3 + 24)--*s 2 * 



and 
~ im JTT 

A-+1 \Ue3 



AP 



U, 



A*3| 



= 2 



+ \U, 



2 m z c 



2 „2 



+ /3 



e3 | 



'2* 



+ 41- 



2^ 4 



4^3 

H — «2tf + W S2* 



,2 



* s 2 * + (5 - 8 ^ 2 ) 4 + s 
(3 + 2 c 2 ) 



,2 
'2* 



3 Can the NuMI-MINOS Combine See This? 

With the Fermilab Main Injector serving as a neutrino source, the MINOS detector [18] at 
the Soudan Mine offers a baseline of 732 km and hence a possible setup for making such a 
measurement. This was investigated by Choubey and Roy in Ref. [15]. Working with a 
exposure of 16 x 10 20 primary protons on target, they conclude that a significant rate asymmetry 
(AN) can be observed for moderate values of (1 — 2\U fJi3 \ 2 ) and |[/" e3 | 2 . 

However, while their idea was an interesting one and the formalism essentially correct, a 
reestimation of the rate asymmetry using the same set of inputs reveals some discrepancies. 
While we agree with Ref. [15] on the magnitude of AiV for all the different cases of inputs 
parameters, we completely disagree on the issue of the error bars. Even on using an optimistic 
indicator of the background fluctuation, namely ^N^ + iV^, we find that, for the abovemen- 
tioned luminosity, the statistical error bars are large enough for each of the rate asymetries in 
Ref. [15] to be consistent with zero. 

Thus, to obtain a result of any statistical significance would need the flux (or, equivalently, 
the detector size) to be increased manifold (see Fig. 1). Such a large scaling, unfortunately, 
is not feasible with the present experimental facility. Note, in addition, that this projection 
altogether neglects the systematic uncertainties, and consequently is already too optimistic. In 
reality, several uncertainties are involved in the translation between the visible energy observed 
at a detector and the energy of the neutrino to which the event is ascribed. These include 
the uncertainties in the assumed final state multiplicities, the scattering or absorption of the 
secondary particles etc. In a recent paper [19], the MINERz/A collaboration, estimates the size 
of the corresponding systematic effects at MlNOS-like experiments could be comparable to, or 
even dominate, the statistical errors. 
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Figure 1: 27ie minimal ratio by which the MINOS exposure (16 x 10 20 primary protons on 
target) needs to be scaled up to obtain la effects for each individual energy bin. The curves 
correspond to the four parameter points of Ref. [15]. Note that only the statistical fluctuations 
have been included in this analysis. 



4 The case for a very long baseline 



Having established that, contrary to the claims of Ref. [15], MINOS, in its current incarnation, 
is not sensitive to possible deviations of {U^l from maximality even for optimistic choices 
of parameters, it becomes interesting to investigate if any of the other proposed long-baseline 
experiments would do the job. And, in a similar vein, if even the luminosity scaling as described 
in the previous section would prove adequate. We start by investigating the second problem. 

As has already been emphasised above, it is not only AP that is important, but also the 
average survival probability P av = (P w + P^)/2 for the latter controls the total number of 
events and hence the size of the statistical error (significance oc AP flfl /^2P !W ). In Fig.2, we 
present both these quantities as function of neutrino energy for a baseline of 732 km. The two 
choices of {U^\ 2 , \U e ^\ 2 ) correspond to the most (and least) optimistic sets of parameters used 
in Ref. [15]. 
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Figure 2: (a) The neutrino survival probability as a function of neutrino energy for a baseline 
L = 732 km. The curves correspond to different combinations of Sm 2 3 (in eV 2 ), \U^\ 2 and 
| f/ e 3 1 2 ■ (b) The corresponding difference in the z/ M and survival probablities. 
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It is abundantly clear that, for such a baseline, and for realistic choice of parameters, AP is 
very small, especially at larger neutrino energies. This was to be expected as a baseline of 732 
km implies that the neutrino beam neither travels through a very dense matter core nor does 
it travel through a "dilute" segment long enough for the matter effect to build up sufficiently 
Clearly, it is dangerous to claim a signal based on AP £ 10~ 4 , for such an act presupposes a 
very accurate knowledge of both the (anti-)neutrino cross sections as well as the density profile 
within the earth. Thus, one would have to concentrate only on neutrinos with relatively small 
energies (E u £ 2 GeV) . Apart from the fact that the asymmetry is not too large even in this 
range 7 , such a measurement also necessitates both a very good energy resolution as well as 
a low energy threshold. Furthermore, for such small energies, neutrino cross-section scales as 
E"l and thus restricting ourselves to a small window is tantamount to rejecting a very large 
fraction of the events. 




E v (GeV) E v (GeV) 




The situation, vis. a vis. the size of AP, seems to improve dramatically with an increase 
in the baseline. For definiteness, we consider two particular values, namely L = 4865 km (the 
distance between the Japan Hadron Facility and the proposed Indian Neutrino Observatory 

7 In fact, it stands to be swamped by the systematic errors unless further experiments are performed [19]. 
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site at Rammam [14]) and L = 10480 km (the distance between Fermilab and INO). Since the 
matter term can no longer be treated as a perturbation, we desist from using the simplified form 
of eq.(9) and use the full expression for neutrino propagation in matter with the density profile 
being given by the Preliminary Reference Earth Model [16]. As Fig. 3 amply demonstrates, for 
each of these values, AP is large over a wide range of E v . In fact, for the larger of the two 
baselines |AP| can be as large as 0.8 for certain E u bins. Furthermore, unlike in the case of 
the 732 km baseline where AP was nearly monotonic in E v , significant structure, including 
sign reversal, is shown. As can be easily appreciated, the latter feature has the potential of 
serving as a key experimental signature. 

A further point of interest relates to the position of the maxima of the survival probabilities 
as well as the asymmetry. Clearly, a high resolution measurement of the same could, in 
principle, be used to determine 5m\ z . 

Before we interpret Fig. 3 to imply that the asymmetry measurement in very long baseline 
experiments would be a sensitive probe of the deviation of \U^\ from maximality, it should 
be realised that the latter is not the only source for a non-zero asymmetry. In fact, as has 
already been hinted at above, the higher order terms in A that were dropped while deriving 
eq.(9) are not proportional to (1 — 2|[/ At3 | 2 ). This can also be seen explicitly from eqs.(14 & 
16) obtained assuming a constant density profile for earth. The numerical importance of such 
contributions naturally increases with the baseline, and has been explicitly displayed in Fig. 4. 
In fact, for the larger baseline, it is quite apparent that the bulk of the effect is due to such 
"higher-order" terms. Thus, it may be concluded that while a baseline of ~ 5000 km may 
allow for a determination of (1 — 2\U^\ 2 ) through measurements of the rate asymmetry, for 
baselines much longer than this the sensitivity reduces quite sharply, and a measurement of 
the maximality of \U^\ is not very straightforward for ultra-long baselines. 
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Figure 4: The normalized probability asymmetry for a baseline of (a) 4865 km and (b) 104 80 
km. In each case, the solid (dashed) line refers to [U^ 2 = 0.36 (0.50). Note that, according 
to the leading expression of eq.(9), the dashed curves should have coincided with the abscissa. 

On the other hand, the asymmetry must remain proportional to \U e 2,\ 2 even on the inclusion 
of A n terms, for it is only through the interaction with the v e that matter effect truly enters 
v^-Vr oscillations. Thus, a baseline of 10480 km may still serve the purpose of affording a 
good measurement of \U e ^\. In fact, since and U e s are inextricably linked even for the case 
of the 4865 km baseline, such an independent measurement of the latter element has its own 
importance. 
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4.1 v e — > ^ oscillations 

With the matter effect playing such an important role in the survival probability for the muon 
(anti-)neutrino, it becomes interesting to consider its effect on P[y e —* v^). For, a source that 
produces [p^) copiously, whether it be a beam dump or a muon storage ring, would also 
produce electron (anti-)neutrinos. 
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Figure 5: The energy dependence of the transition probability for v e — > (left) and v e — > P M 
(right) on passage through matter. The upper and lower panels correspond to baselines of 4865 
and 104 80 km respectively. 

In Fig. 5, we present the probabilities for such transitions as a function of (anti-)neutrino 
energy. Several observations are in order: 

• The transition probability P{y e v^) can be quite sizable for parameter values that are 
still allowed by experimental data. 

• For E £ 3 GeV, we see that P{p e — > Vp) ~ P{y e v^) (except, of course, near the node 
at ~ 8 GeV for the 10480 km case). Coupled with the smaller detection efficiency for 
anti- neutrinos, this indicates that it would be more profitable to work with v e than u e . 



The transition probabilities have non-trivial dependence on both U e 3 and U^. Unlike 
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in the case of the aforementioned asymmetry, these are not in general proportional to 

\u e3 \ 2 . 

• The measurement of this effect would thus lead to a constraint in the U e3 -U^ plane 
independent of that drawn from the asymmetry measurement. 



5 Event rate calculation from a storage ring neutrino 
source 

As we have already mentioned in the introductory section, the neutrino flux from a muon 
storage ring can be calculated very precisely. Starting with a pT beam, the number of /i~ 
events in a far detector can be obtained by folding this flux with the survival probability 
P MM (L, E) and the charged current cross-section: 

r dN 
N, = N n J <j(iv + N - y~ + X) P^(L, E) dE v ^ (17) 

where N n is the total number of nucleons present in the fixed target. An analogous expression 
obtains for /i + events as well. Note though that the v e from muon decay could also oscillate 
into while traversing through the earth and result in muonic events. However, they result 
in wrong sign muons and thus can be easily distinguished in a magnetized detector. We shall 
revert back to them at a later stage. 

To be specific, we shall consider a storage ring with 20 GeV muons and a 50 kT Iron 
calorimeter detector one such as the proposed monolith [13] or the ical/ino experiment [14]. 
The projected energy threshold for muon detection is about 2 GeV and the resolution is 
expected to be better than 0.5 GeV over the entire range. 

Before delving into the actual event profile, several points are in order. 

• As has already been demonstrated in the previous section, the difference in the survival 
probalities (AP MM ) typically grows with the baseline. Thus the experimental sensitivity 
could be expected to increase with the baseline. (This is modulo the fact that the 
sensitivity of the asymmetry to (1 — 2 1 1 2 ) reduces for very large baselines, while the 
sensitivity to \U e3 \ 2 is retained.) 

• On the other hand, for a fixed detector size, the effective solid angle subtended by the 
detector to the source decreases with increasing baseline. Consequently the neutrino 
flux goes down quadratically with the baseline, thereby resulting in a smaller number of 
events and hence larger statistical errors. This, for example, can be read off from Fig. 6, 
wherein we have displayed the number of events expected if the survival probability were 
to be unity. 

• One obvious way to circumvent the latter problem is to start with higher luminosity. For 
our numerical results we then adopt the 'higher luminosity' option for the muon beam, 
namely 10 20 muon decays per year. 

• An alternative way to the same end would have been to start with a higher energy for the 
muon beam. For example, the neutrino beam from a 50 GeV storage ring is collimated 
sufficiently enough to get a similar significance with even the 'low luminosity' option for 
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Figure 6: The number of events expected as a function of the (anti-) neutrino energy if the 
survival probability were to be unity. The source is a 20 GeV muon storage ring with 10 20 
(anti-)muons decaying while the detector is a 50 kT iron calorimeter with a energy threshold 
of 2 GeV for the ^ [14]- The left (right) panel corresponds to a baseline of 4865 (10480) 
km. Also shown are the events generated by the electron (anti-) neutrino assuming a transition 
probability of unity. 



the storage ring (10 19 muon decays per year). However, we shall desist from using this 
option as such a machine is likely only as a second generation facility. 

• As Fig. 3 amply demonstrates, for longer baselines such as the one we are considering, 
the survival probabilities are very sensitive functions of the (anti-)neutrino energy, espe- 
cially for low energies. Thus, extracting any information from the low energy tail of the 
detector muon spectrum would necessitate very good energy resolution. However, with 
the neutrino- nucleon cross sections being small for such E u , these neutrinos have only a 
relatively small contribution to make to the total number of muon events. 

Keeping in mind that v — N cross-section is nearly double 8 that of v — N we argue that the 
detector should be exposed to a beam for double the time that it is exposed to a beam. 
To be specific, we consider a 50 kT-year exposure to a z/ M beam and a 100 kT-year exposure 
to beam, or in other words, a total exposure of 3 years for the detector configuration [14] 
under consideration. 



5.1 The rate asymmetry for v^iy^) initiated events 

Armed with the above, we can now calculate the difference in \i~ and fi + event numbers. 
Defining an asymmetry of the form 

N- - N, 



A N = 



8 For the total cross section, this ratio has only a mild dependence on the incident neutrino energy, and 
for the energy range we are interested in has a value close to 2.18. However, on imposition of an energy 
threshold £?(/x ± ) > 2 GeV that such a detector is expected to have, the ratio is a more sensitive function of 
the (anti-)neutrino energy (see Fig. 6). In our numerical analysis, we explicitly account for this factor and the 
consequent effect on the error bars. 
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we present, in Fig. 7, the same for a baseline of 4865 km. The values for (jEV 2 , | f^e3 1 2 ) 
correspond to the two extreme cases of Ref. [15]. Also plotted are the analogous expectations 
for \U^\ 2 = 0.5, a case that, to the lowest order, would have been associated with vanishing 
An. We have binned the expected data, keeping under consideration both the expected energy 
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Figure 7: The difference between and charged- current events for a baseline L = 4865 
km with an integrated luminosity of 10 20 muons and double the number of antimuons. The 
detector is a 50 kT iron calorimeter with a muon energy threshold of 2 GeV [14]- Each plot 
corresponds to a different combination of (Sm 2 3 , |EV 2 ; lE/^! 2 ,). The other parameters are: 



Sm 2 12 = 7 x 10~ 5 eV 2 and sin 2 12 



0.3. Also shown are the la error bars. 



resolution [14] as well as the number of events in a particular bin. Also shown are the statistical 
la error bars. Several features are easily discernible: 

• In each of the cases displayed, the rate asymmetry is statistically significant in more 
than one bin. Together, they lead to a very considerable shift in the x 2 ■ This amply 
demonstrates the sensitivity of such an experiment to a departure from maximality of 

• As expected, the signal is less pronounced for the case (|EV 2 , |E/" e 3| 2 ) = (0.40, 0.05) than 
for (0.36, 0.07). This is just a reflection of the fact — see eq.(9) — that the rate asymmetry 
is proportional to |E/ e 3| 2 . While the proportionality to (1 — 2 |EV 2 ) is no longer an 
accurate statement, it, nevertheless still encapsulates a large measure of truth. 
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• In the energy range of interest, the binwise asymmetry clearly shows an oscillatory be- 
haviour, a feature that would have been absent for a similar detector had the baseline 
been shorter than ~ 1000 km. 

• The functional dependence of AN(E U ) exhibits a discernible dependence on the mixing 
angles (|?7 M 3| 2 , |t/ e3 | 2 ), and could, in principle, be used to distinguish between values for 
these parameters. However, to obtain a significant resolution in the parameter space, a 
much larger event count would be necessary. 

• The dependence on the value of 5m\ z , on the other hand, is much more pronounced. 
The shift in AP that we encountered in Fig. 3 is quite well reflected by a shift in An, even 
after the convolution with the muon spectrum and the energy-dependent cross sections 
as well as the finite resolution effects. Thus, our previous comment about using this 
measurement for a determination of 5m 2 3 is substantiated. 

• The energy dependence is quite different for the two baselines considered and can be used 
to advantage. For example, if storage rings come up at both the JHF [20] and Fermilab, a 
detector such as ical/ino [14] could use beams from both to distinguish more efficiently 
between the possible parameter sets. 

5.2 Wrong sign muons 

We now turn to a discussion v e to conversion (as well the conjugate process) and their 
detection in the iron calorimeter detector under consideration. Such events are obviously 
characterised by the appearance of a muon with a charge opposite to that of the decaying 
particle in the storage ring. Since the detector is to be magnetized one, charge measurement 
is relatively straightforward thus rendering these events easily distinguishable. Efficacy of 
'wrong-sign-muon' events in measuring the angle #13 has already been advocated in Ref. [10] 
in the context of muon-storage ring neutrino factory experiments. Here we want to emphasise 
that how an independent measurement of 'wrong-sign-muon' events can also supplement the 
more accurate measuremeant of # 2 3- We have already seen, vide Fig. 5, that the conversion 
probability for u e is typically larger than that for v e , a consequence of the sign of 5m\ 2 - This 
fact, alongwith the smaller detection efficiency for v^, implies that it is wiser to concentrate on 
the fi + — > v e — > chain rather than attempt to measure the conjugate process as well. And 
since we have already advocated twice as long an exposure to fi + decays as compared to pT 
decays, this works in our favour too. 

In Fig. 8, we present the number of \i~ events in the detector when exposed to ji + decays. 
As before, the detector is supposed to be a 50 kT one, with a total exposure of 2 x 10 20 decays. 
Clearly, the number of events grows with |£/ e 3| 2 . While the exact dependence is not linear (as 
the naive expectation would be), explicit calculations with other values of \U e ^\ shows that it 
is not very far from being linear. This is also the case for the larger baseline of 10480 km 
(which we do not exhibit here). The dependence on \U^\ is more complicated though. As 
can be expected, this depends more crucially on the size of the matter effect and hence on 
the exact baseline (see, for example, Fig. 5). However, given the considerably large number of 
such events, it is quite apparent that this measurement is likely to be a sensitive probe in the 
( I j , \U e z\) plane. Unfortunately though, the event distribution is not very sensitive to <5m 2 3 . 

The main motivation for the above excercise is to show that the measurement of any of the 
small parameters (\U e s\ or j3) at long baseline experiments cannot be done independent of each 
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Figure 8: The number of pT events in a 50 kT iron calorimeter exposed to a total of 2 x 10 20 p + 
decays in a 20 GeV storage ring at a distance of 4865 km. The histograms correspond to differ- 
ent values of (It^l 2 , |£4s| 2 ) with the left (right) panels referring to 5ml 3 = 0.002 (0.0027) eV 2 . 



other. Although experiments with relatively small baselines (such as NuMi-MlNOS) seem to be 
quite sensitive to the parameter (3, this optimism is misplaced as explained in Sections 3 & 4. 
On the other hand, moderate baseline experiments do show increased sensitivity to (3, but only 
at the cost of having a more complicated dependence on (3 and |C/ e3 | (see eq.14 and Fig. 4). 
In either case, a measurement of (3 from ^-survival probability asymmetry needs an accurate 
knowledge of \U e s\. It is here that the wrong-sign muon rates (which is more sensitive to U e s) 
have the most important role to play. Of course, a precise knowledge of Smf 3 is important 
for the extraction of \U e3 \ (although the dependence on 8m\ 3 is not as pronounced for the 
wrong-sign muon events as for the right-sign ones). However, it is expected that the error-bars 
on 5m\ 3 would be significantly reduced before the era of such neutrino factory experiments. 
Then, in principle, An and wrong-sign muon rates could be used simultaneously in a statistical 
analysis (in a spirit similar to that followed in Ref. [11]) to extract \U e ^\ and (3 to much greater 
accuracy than possible today. 



6 Conclusions 

Analyses of matter effects in oscillations have largely concentrated on transition proba- 
bilities. Revisting the problem, we demonstrate that the survival probabilities too are very 
sensitive to matter effects and can be used profitably in determination of several crucial pa- 
rameters in the neutrino sector. 

To start with, we derive a set of approximate analytical expressions for the survival 
probability in presence an arbitrarily large matter density. While similar results had already 
been reported for small matter terms in the effective Hamiltonian, our expressions have a much 
wider range of validity. 

Since we are interested in the survival probability for v^, a variable of interest can be 
constructed simply from the number of charged current (CC) events in the detector. Starting 
with and beams, A N = (iV_ — N + )/(N_ + iV + ) describes an asymmetry between the 
number of 'same-sign' p~ and p + events generated in the detector after traversing a given 
baseline. We show by explicit event rate calculation that this asymmetry is a good measure of 
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the matter effects felt by the neutrinos while propagating through the earth. However, contrary 
to claims in the literature [15], the size of the asymmetry for the Fermilab-MINOS combine is 
much smaller than the experimental sensitivity. 

We suggest, therefore, that a future long-baseline experiment could explore this effect. As 
a prototype experiment, we consider a neutrino factory (a muon storage with muon beams of 
energy 20 GeV) and the proposed 50 kT iron calorimeter detector [13, 14] with a capability 
of muon charge determination. Using realistic experimental resolutions, we demonstrate that 
such a combination is sufficient to establish the aforementioned asymmetry with a very large 
statistical significance. However, unlike in the case of relatively short-baseline Fermilab-MINOS 
combine, the asymmetry for a long baseline is no longer proportional to the deviation of \U^\ 2 
from maximality. Even so, the sensitivity to (1 — 2 |t/^3| 2 ) remains quite pronounced for 
baselines of upto about 6500 Km. Using the aforementioned approximate expressions for the 
survival probabilities, such data can thus be used to determine (1 — 2|[/ M 3| 2 ). While the 
sensitivity to this combination decreases for much longer baselines, asymmetry data for such 
baselines can still be used for precise determination of \U e3 \. 

The latter exercise is also shown to be aided by the measurement of 'wrong-sign-muon' 
events generated through v e — > z/ M oscillations. Estimating the rate of such events for the same 
experimental setup, we demonstrate that this could be used in conjunction with the asymmetry 
measurement to lead to a good simultaneous determination of \Ufjz\ and \U e ^\. 
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